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Abstract — Simple algorithms for retrieving free-space antenna 
field or directivity patterns from complex (field) or real (intensity) 
measurements taken in ideal reverberation environments are 
introduced and discussed. 

I. Introduction 

Antenna measurements are usually performed under simu- 
lated free-space conditions, e.g. by placing the antenna under 
test (henceforth AUT) as well as the measuring probe in 
an open test-range or in an electromagnetic (henceforth EM) 
anechoic chamber [1]. 

Retrieving free-space antenna field and/or directivity pat- 
terns from measurements taken in any realistic (i.e., imperfect) 
test-range or anechoic chamber relies on the possibility of re- 
constructing the ray skeleton of the measured field using robust 
spectral estimation techniques, including, e.g., periodogram, 
Prony, Pisarenko, Matrix-Pencil and Gabor algorithms [2]-[5], 
so as to "subtract" all environment-related reflected/diffracted 
fields. 

The direct (free-space) field, however, cannot be unambigu- 
ously identified, unless additional assumptions are made about 
its relative intensity and/or phase, which do not hold true in 
the most general case. 

A possible way to uniquely extract the free-space (direct 
path) field is to average over many measurements obtained by 
suitably changing the position of the source-probe pair with 
respect to the environment, while keeping the source-probe 
mutual distance and orientation fixed. This obviously leaves 
the free-space direct-path term unchanged, while affecting 
both the amplitudes and the phases of all environment-related 
reflected/diffracted fields. In the limit of a large number 
of measurements, one might expect that these latter would 
eventually average to zero. This is the rationale behind the 
idea of retrieving free-space antenna parameters from mea- 
surements taken in a reverberation enclosure (henceforth RE), 
where the chamber boundary is effectively moved through 
several positions by mechanical stirring, while the source- 
probe distance and mutual orientation is fixed. 

Through the last decades reverberation enclosures earned 
the status of elicited EMI-EMC assessment tools [6]. On the 
other hand, only recently effective procedures for estimating 
antenna parameters, including efficiency [7], diversity-gain [8], 
MIMO-array channel capacity [9], and free-space radiation 
resistance [10], from measurements made in a reverberation 
chamber have been introduced by Kildal and co-workers, in a 
series of pioneering papers. 
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Here we discuss, perhaps for the first time, free-space 
antenna field/directivity pattern retrieval from measurements 
taken in a reverberation environment. 

The paper is organized as follows. In Sect. |ll] the key 
relevant properties of reverberation enclosure fields are sum- 
marized. In Sect. [Hi] and IIVI simple algorithms for retrieving 
free-space antenna field or directivity patterns, respectively 
from (complex) field or (real) intensity measurements made 
in a reverberation environment are discussed, including the 
related absolute and relative errors. The related efficiency is the 
subject of Sect. [V] including some useful concepts on Cramer- 
Rao bounds. Conclusions follow under Sect. IVII 

II. Fields in Reverberation Environments 

In the following we shall sketch and evaluate some straight- 
forward procedures to estimate the free-space antenna field or 
directivity pattern from measurements made in a reverberation 
enclosure. 

The AUT field/intensity will be sampled at a suitable 
number of points P of the AUT-centered sphere r = R, 
corresponding to as many sampling directions. At each point 
P we shall actually make N measurements in the reverberation 
environment corresponding to as many different positions of 
the mode stirrers. 

Throughout the rest of this paper we shall restrict to the 
simplest case where both the antenna under test and the field- 
probe (henceforth FP) are linearly (co)polarized, and placed 
in an ideal (fully-stirred) reverberation environment. 

The relevant component of the complex electromagnetic 
(henceforth EM) field at a point P can be written: 

E(P,n)=E d (P)+E r (P,n), n=l,2,...,N. (1) 

The first term in Q is the direct field, and is the only 
term which would exist in free-space; the second term is the 
(pure) reverberation field, whose value depends on the stirrers' 
positions 1 , and n identifies the different stirrers' positions. 

According to a widely accepted model [11], for any fixed 
P, the set {E r (P, n)\n = 1,2,..., N}, can be regarded as an 
ensemble of identically distributed (pseudo) random variables 
resulting from the superposition of a large number of plane 
waves with uniformly distributed phases and arrival directions. 

Under these ideal (but not unrealistic) assumptions the real 
and imaginary part of the reverberation field E r (P, n) will be 
gaussian distributed 2 and uncorrelated, with zero averages and 

'We consistently include in the free-space antenna-field any re- 
flected/diffracted term which does not change as the positions of the mode 
stirrers change over. 

2 In this connection, the amplitude distribution of the contributing plane 
waves turns out to be almost irrelevant [12]. 
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equal variances [13], 



(Re 2 E r (P,n)) = {lm 2 E r (P,n)} 



E 2 



(2) 



where (•) denotes, more or less obviously, statistical averaging. 
The quantity E 2 in <|3J is given by [11]: 



8nm 

A 2 



(3) 



where r/o is the free-space wave impedance, A the wave- 
length, and n r the power received by any (linearly polarized, 
matched) antenna placed in the reverberation enclosure, irre- 
spective of its orientation and directivity diagram [11]. This 
latter is related to the total power n t radiated into the enclosure 
by the AUT as follows [14] : 



(4) 



where the (frequency dependent) RE calibration-parameter 3 
is related to the chamber (internal) surface S and wavelength 
A by [14] 

A 2 

a being an average-equivalent wall absorption coefficient 3 . 



III. AUT Free-Space Field Estimator 

Under the made assumption where the real and imaginary 
part of the reverberation field E r are independent, zero-average 
gaussian random variables, it is natural to adopt the following 
estimator of the free-space (complex) AUT field at P in terms 
of the (complex) fields Q: 



E d (P) = N- 



N 

£ 

71=1 



E{P, n) 



(6) 



Equation (|6jl provides unbiased estimators of Re[E d (P)] and 
lm[Ed(P)], with variances 



VAR[ReE d {P)} = VAR[lmE d (P)} = ||. 



(7) 



The related absolute and relative errors are: 



and 



e (F) 
^abs 



-rel 



E d {P) - E d (P) 



E d {P) - E d (P) 



1/2 



distance fixed, this will at the same time make 9 larger, in 
view of eq. dlOi . thus reducing the relative error (|9j, when 
meaningful, as well. Note that this is true for both far and 
near-field measurements. 

IV. AUT Directivity Estimator 

The AUT directivity can be estimated from (far field) 
intensity measurements made in a reverberation enclosure as 
follows. Let 

I(P,n)=Re 2 E(P,n)+lm 2 E(P,n) : n—1, 2,. . ., N. (11) 

It is convenient to scale the field intensities I(P, n) to the 
variance in 0, by letting £„(P) = 2I(P,n)/E 2 , so that all 
the £„ are (identically) distributed according to a noncentral 
chi-square with two degrees of freedom [15] and non-centrality 
parameter 9{P) given by eq. dlOi . 

We may use the obvious far field formula: 



i^)i 2 = g§ D in 



(12) 



where D(P) is the AUT directivity, together with eq.s Q and 
10} in eq. (tlPl to relate to the AUT directivity D as follows 



1 



A 



D =: jD, 



(13) 



8tt 2 5(/) \R, 

where the dependence of D and 8 on the measurement point 
(direction) is understood and dropped for notational ease 4 . 
The probability density function of the can be accordingly 
written [15] 

m D) = 1 exp (-^y^-) ib(VfrD). (14) 
The first two moments accordingly are: 

(£)=2 + 7.D, VAR{€]=4: + 4rfD, (15) 

which suggest using the following (simplest, unbiased) esti- 
mator of D [16]: 



= A^ 1/2 E Q , (8) f or W hich 



D = (7 AT) 



VAR[D] 



N 

_1 £ 

71 = 1 



(&. - 2), 



j 2 N 



(1+7-D). 



(16) 



(17) 



1/2 



2 \ 1/2 
= 'M ' (9) 



The absolute and relative errors of the directivity estimator 
H6i are thus: 



\E d (P)\ 

where, for later convenience, we introduced the dimensionless 
quantity 



e (D) 
abs 



VAR\D] 



and: 



E 2 

The r.m.s. absolute error l|8} can be made as small as one 
wishes, in principle, by increasing N, and/or the chamber size 
(the distance between the chamber walls and the AUT-FP pair), 
which makes Eq oc S -1 / 2 smaller. Keeping the AUT-FP pair 

3 The coefficient a in |5j can be evaluated as a = Eo/(E + En) < 1, So 
being the area of a wall aperture which halves n r [14]. 



N 



1/2 



(18) 



(19) 



The absolute and relative (when meaningful) errors ( I18> and 
( II 91 can be made as small as one wishes, in principle, by 
increasing N, and/or the chamber size, so as to make 9 suitably 
large, in view of dl 3i . 

4 For the simplest case of a spherical enclosure of radius R c , from eq.s J5} 
and (Uj one gets 9 = (a/ir)(R c / R) 2 D, R being the AUT-FP distance. 
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Note that in all derivations above we made the implicit 
assumption of dealing with independent measurements. 

The number N of independent measurements needed to 
achieve relative errors ~ 5 • 1CP 2 for both field and directivity 
measurements is shown in Fig. 1, and is of the order of 10 2 




100. 



Fig. 1. Number TV of independent measurements vs. 9 needed for relative 
errors ~ 5 ■ 10 -2 . Solid line: e dashed line: e , . 



for 9 ~ 10. This figure is consistent with typical experimental 
findings [17], and also with theoretical estimates obtained from 
a chaos-based models of reverberation enclosures [18]. 

V. Efficiency of Proposed Estimators 

An obvious question is now whether one could do better 
using different estimators, other than (|6j and (I16> . 

The natural benchmark for gauging the goodness of an esti- 
mator is the well-known Cramer-Rao lower bound (henceforth 
CRLB) [19]. We limit ourselves here to remind a few basic 
definitions and properties. Let £ = {£,k\k = 1,2,..., K} a 
set of (real) random variables with joint probability density 
/(&X), where X = {X m \m = 1,2,..., A/} is a set of 
(unknown, real) parameters to be estimated. One can prove 
that 5 for any estimator X of X such that (X) = X, (unbiased 
estimator), one has: 



C x -J~ 1 (X)>0 
where C is the covariance matrix, viz.: 



hk- 



[J(X)] 



hk 



-X h )(Xk-X k ) 

G> 2 log/(£;X) 
dX h dX k 



(20) 



(21) 



(22) 



is the Fisher information matrix, the expectations are taken 
with respect to /(£;X), and the true value of X is used 
for evaluating (I22> . Equation i20\ implies M inequalities 
whereby the diagonal elements of C^, i.e., the variances of 

5 We implicitly assume that the following regularity condition [15] holds: 



the components of X, are bounded from below. These are 
the CRLB s. An estimator for which the l.h.s. of eq. ( I20t is 
actually zero, i.e., for which the variance of each component 
of X attains its CRLB is called efficient. For the special case 
where the £/. are independent and identically distributed, with 
a PDF /(£; X) depending on a single parameter X, equation 
(I20> becomes 



VAR[X] > -- 



1 



d 2 log/(£;A) 



N 



dX 2 



N i 



dlogmx) 



dX 



(23) 



One can readily prove that the field estimator (|6} is an 
efficient one, since the r.h.s of @ coincides with the pertinent 
CRLB. The simplest directivity estimator dl6> . on the other 
hand, while not efficient, gets very close to its CRLB, as shown 
below. 

The Cramer-Rao bound for the estimator d!6i . is obtained 
by using the following formula, which follows directly from 
eq. (HI, 



d log /(£;£>) 



2 

2 



1 h (V&D) 



3D 2 2 1 (Vf^D) 

where 2o,i(") are modified Bessel functions, and is [16]: 

4 / A 
CRLB(D) = _ _ - 1 



j 2 N \jD 



(24) 



(25) 



where 



A = 



In 



(26) 



and the expectation is taken with respect to /(£; D). 

The ratio between the CRLB J26I and the variance Mil 
yields the relative efficiency 



CRLB(D) 



(27) 



VAR[L>] 

The relative efficiency i27\ of the proposed directivity esti- 
mator dl6> is readily computed from eq.s i25\ . ( I26t and d!7l >. 
and is independent of N. It is displayed in Fig. 2 vs. 8 = jD. 
The relative efficiency of d!6l > is seen to be pretty decent, being 
always larger than « .937. 

VI. Conclusions 

Free-space antenna field/directivity measurements in ideal 
reverberation enclosures have been shortly described and 
evaluated. The main simplifying assumptions (linearly co- 
polarized AUT and FP) can be more or less easily relaxed at 
the expense of minor formal complications which do not alter 
the main conclusions. On the basis of these preliminary results, 
the possibility of performing cheap, simple and reliable in situ 
antenna measurements using, e.g., flexible conductive thin-film 
deploy able/inflatable enclosures with air-blow stirring [20], 
[21] is envisaged. 
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Fig. 2. Directivity estimator. Relative efficiency p, eq. 1271 . vs. 9, eq. 1131 . 
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